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Abstract 



On the basis of the nonlinear theory of elasticity, the general constitutive equa- 
tion for an isotropic hyperelastic solid in the presence of initial stress is derived. 
This derivation involves invariants that couple the deformation with the initial 
stress and in general, for a compressible material, it requires 10 invariants, reduc- 
ing to 9 for an incompressible material. Expressions for the Cauchy and nominal 
stress tensors in a finitely deformed configuration are given along with the elasticity 
tensor and its specialization to the initially stressed undeformed configuration. The 
equations governing infinitesimal motions superimposed on a finite deformation are 
' then used to study the combined effects of initial stress and finite deformation on 

s * . the propagation of homogeneous plane waves in a homogeneously deformed and 

initially stressed solid of infinite extent. This general framework allows for vari- 
ous different specializations, which make contact with earlier works. In particular, 
■ connections with results derived within Biot's classical theory are highlighted. The 

general results are also specialized to the case of a small initial stress and a small 
pre-deformation, i.e. to the evaluation of the acoustoelastic effect. Here the formu- 
las derived for the wave speeds cover the case of a second-order elastic solid without 
initial stress and subject to a uniaxial tension [Hughes and Kelly, Phys. Rev. 92 



(1953) 1145] and are consistent with results for an undeformed solid subject to a 
residual stress [Man and Lu, J. Elasticity 17 (1987) 159]. These formulas provide a 
basis for acoustic evaluation of the second- and third-order elasticity constants and 
of the residual stresses. The results are further illustrated in respect of a prototype 
model of nonlinear elasticity with initial stress, allowing for both finite deformation 
and nonlinear dependence on the initial stress. 
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1 Introduction 



The presence of initial stresses in solid materials can have a substantial effect on their sub- 
sequent response to applied loads that is very different from the corresponding response 
in the absence of initial stresses. In geophysics, for example, the high stress developed 
below the Earth's surface due to gravity has a strong influence on the propagation speed 
of elastic waves, while in soft biological tissues initial (or residual) stresses in artery walls 
ensure that the circumferential stress distribution through the thickness of the artery wall 
is close to uniform at typical physiological blood pressures. Initial stresses may arise, for 
example, from applied loads, as in the case of gravity, processes of growth and develop- 
ment in living tissue or, in the case of engineering components, from the manufacturing 
process, either by design to improve the performance of the component or unintentionally 
(which can lead to defective behaviour). Here we use the term initial stress in its broadest 
sense, irrespective of how the stress develops. This includes situations where the stress is 
due to an applied load leading to an accompanying finite deformation from an unstressed 
configuration, in which case the term prestress is commonly used, and situations in which 
the initial stress arises from some other process, such as manufacturing or growth, and 
is present in the absence of applied loads. If an initial stress is present in the absence of 
applied loads (body forces and surface tractions) it is referred to as residual stress, as in 
the definition adopted by iHoger (jl985h . 

In the context of finite deformation elasticity theory much motivation for the study of 
initial stresses, more especially residual stresses, comes from soft tis sue biomechanics, and 
for a r ecent discussion of residual stresses in artery walls we refer to iHolzapfel and Ogden 
(120101 ) and references therein. The development of constitutive laws for resid ually stressed 



materials h as b een a focus of much of the work of Hoger, as exemplifie d by iHogerl (11985, 
1986 . 1993 ) and Johnson and Hoger ( 1993 ); see al so the recent paper b y Saravanan ( 20081 ). 
In the context of acoustoelasticity the work of Man and Lul ( 19871 ) was based on the 
developments of Hoger, and it also relates to the much earlier work of Biot. Biot's 
work was developed within the geophysical context and he was particularly concerned 
with the effect of initial stress on the propagation of small amplitude elasti c waves. H e 
developed a static theory of small deformations influenced by in i tial st ress flBiotl . Il939h . 



followed by a corresponding theory for wave p ropagation (IBiotl . Il940l ). all of this work 



being conveniently collected in his monograph (IBiotl . Il965[ ) . 



Whatever the source of the initial stress and whether or not it is accompanied by a 
finite deformation, of particular interest in many applications is the effect that the ini- 
tial stress has on small deformations (static or time- dependent), specifically deformations 
linearized relative to the initially-stressed state and often referred to as incremental de- 
formations or motions. If there is an accompanying finite deformation then the theory 
is often referred to as the theory of small deformations superimposed on large deforma- 
tions. This theory requires knowledge of the associated elasticity tensor, which in general 
depends on any finite deformation present and on the initial stress. Thus, part of the 
purpose of the present work is to obtain the general form of the elasticity tensor for an ini- 
tially stressed material that is subject to a finite deformation where the initial stress is not 
itself associated with an initial finite deformation. We shall be concerned primarily with 
the case in which the material possesses no intrinsic anisotropy so that any anisotropy 
arises solely from the presence of initial stress, and we show how the components of the 
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elasticity tensor relate to those of the elasticity coefficients in Biot's isotropic theory. In 
the present paper we use the elasticity tensor to examine the effect of initial stress on the 
propagation of small amplitude elastic waves for both finitely deformed and unstrained 
materials, with reference to the associated strong ellipticity condition. 

In Section [21 the basic equations for a (hyper) elastic material subject to initial stress 
are summarized together with the equations governing equilibrium of the initially stressed 
and finitely deformed initially stressed configurations. The equations of incremental mo- 
tions superimposed on an initially stressed and finitely deformed equilibrium configuration 
are then given in Section [3j wherein the elasticity tensor is derived and its general sym- 
metry properties recorded. In Section 2] we then specialize the general theory and focus 
on the development of the constitutive law for an initially stressed material that has no 
intrinsic anisotropy such as might be associated with preferred directions, for example, 
and any anisotropy in its response relative to the initially stressed configuration is entirely 
due to the initial stress. The constitutive law of the material is based on a strain-energy 
function (defined per unit reference volume) that depends on the combined invariants of 
the right Cauchy-Green deformation tensor and the initial stress tensor. For a compress- 
ible material there are 10 such independent invariants in the general three-dimensional 
case, a number which reduces to 9 for an incompressible material. Expressions for the 
Cauchy stress and nominal stress tensors and the elasticity tensor are given in general 
forms for both compressible and incompressible materials and then specialized for spe- 
cific applications by reducing the number of invariants involved while ensuring that the 
effects of initial stress are adequately accounted for. Some details of the calculations are 
relegated to Appendix A for convenience of reference. 

In Section [5j the equations of motion are specialized in order to study the effect of 
initial stress on the speed of infinitesimal homogeneous plane waves. First, we give results 
for the situation in which there is initial stress but no finite deformation. It is noted, in 
particular, that the wave speed depends in a nonlinear fashion on the initial stress. Re- 
sults are compared with those arising from Biot's isotropic theory based on connections 
between the components of the elasticity tensor used here and those used in Biot's theory. 
The general conne c tion b etween these components, which was derived in a recent paper 
(jOgden and Singhl . l201ll ). is noted here for reference. Also, within the p resent frame- 
work, by considering a linear elastic material, we confirm a formula given by lMan and Lu 
(119871 ) in the general linear theory, which they attribute to Biot, concerning the effect 
of initial stress on the speed of homogeneous plane waves and the acoustoela stic effect. 
Next, for the Murnaghan form of strain-energy function (iMurnaghanl . 119371 ). which is 
appropriate for second-order elastic deformations, a nd specializing to t he ca se of a ma- 
terial without initial stress, we recover the results of iHughes and Kellyl ( 119531 ) concerned 
with the second-order correction for the speeds of longitudinal and transverse waves in an 
isotropic elastic material. The final illustration introduces a prototype nonlinear model of 
an elastic material with initial stress that allows for both finite deformation and nonlinear 
initial stress. 
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2 Basic equations for an elastic solid with initial stress 



Consider an elastic body in some reference configuration, which we denote by B r . Let 
X be the position vector of a material point in B r . Any subsequent deformation of the 
body is measured from B r and we assume that there is an initial (Cauchy) stress r in 
this configuration. This initial stress is symmetric (rotational balance in the absence of 
intrinsic couple stresses) and satisfies the equilibrium equation 

Divr = 0, (2.1) 

in the absence of body forces, where Div denotes the divergence operator with respect 
to B r . If the traction on the boundary dB r of B r vanish es pointwise t hen r i s refe rred 



to as a residual stress, and it is necessarily non-uniform flHogerl . Il985l : lOgdenl . 120031 ). If 
the traction is not zero then the initial stress may or may not be accompanied by some 
prior deformation required to reach the configuration B r from a completely unstressed 
configuration. Here we shall not be concerned with how the initial stress is produced. 

Suppose now that the body is deformed quasi-statically into a new configuration B 
with boundary OB so that the material point X takes up the new position x given by 
x = %(X), where the vector function \ defines the deformation for X 6 B r . The so- 
called deformation x is required to be a bijection and to possess appropriate regularity 
properties, which need not be made explicit here. The deformation gradient tensor, 
denoted F, is defined by F = Grad%, where Grad is the gradient operator with respect 
to B r , and the left and right Cauchy-Green deformation tensors are defined by 

B = FF T , C = F T F, (2.2) 

respectively. 

Let a and S denote the Cauchy stress tensor and the nominal stress tensor, respec- 
tively, in the configuration B. For equilibrium in the absence of body forces a and S 
satisfy the equations 

div<r = 0, DivS = 0, (2.3) 

and we note the standard connection a = (detF) _1 FS. For a material without couple 
stresses, a is symmetric and hence we have 

FS = S T F T . (2.4) 

For an elastic material we consider a strain-energy function W defined per unit volume 
in B r . This function depends on the deformation gradient F and the initial stress r. If r 
depends on X, as would be the case for a residually-stressed material, then the material 
is necessarily inhomogeneous, but if r is independent of X the material is homogeneous 
unless its properties depend separately on X. In either case we make the dependence on 
r explicit and write 

W = W(F,t). (2.5) 

Of course, by objectivity, W depends on F only through C = F T F, but otherwise this 
form of W is completely general and no material symmetry is invoked. Note, however, 
that in general, i.e. if r is not a hydrostatic compression or tension, the presence of r 
induces some anisotropy in the material, even if the material has no intrinsic anisotropy. 
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Thus, r has an effect on the constitutive law analogous to that of a structure tensor in 
anisotropic elasticity. 

For a material not subject to any internal constraints, the nominal and Cauchy stresses 
are given by 

S = ^(F, r), cr = J-'FS = ^F^(F, r), (2.6) 
respectively, where J = detF > 0. When evaluated in B r , these give the connection 

dW 

-=OTrd.-). (2-7) 

where I is the identity tensor. 

For an incompressible material, the internal constraint 

J = detF = l (2.8) 

must hold for all deformations and the counterpart of (I2.6P in this case is 

dW , dW 

S = —(F,T)-pF~\ <x = FS = F— (F,r)-pl, (2.9) 

where p is a Lagrange multiplier associated with the constraint. When evaluated in the 
residually stressed configuration, these both reduce to 



<9F 

where p^ r > is the value of p in B r . 



1,t (2.10) 



3 Incremental motion superimposed on an initially 
stressed configuration subject to finite deforma- 
tion 

Superimposed on the equilibrium configuration B defined by x = x(X), we now consider 
an incremental motion x(X, t), where t is time. Here and in the following a superposed 
dot indicates an incremental quantity, increments are consider 'small', and the resulting 
incremental equations are linearized in the increments. Thus, x represents the displace- 
ment from x, and we shall also express it in Eulerian form by writing the displacement 
vector as a function of x and t, namely u = u(x, t). The corresponding increment in the 
deformation gradient, F, is expressible as 

F = LF, (3.1) 

where L = gradu is the displacement gradient. 

The (linearized) incremental nominal stress takes the forms 

S = <! " 4F ' (3 2) 

' AF + pF^FF- 1 - pF' 1 , 
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for unconstrained and incompressible materials, respectively, where 



d 2 W d 2 W 
A = WdF' Amfij = W~dF^' (3 ' 3) 

is the elasticity tensor and, in component form, AF = AaipjFjp defines the product 
used in the above. Here and henceforth the standard summation convention for repeated 
indices applies. We note in passing that ( 13. 2 p is applicable even if a strain-energy function 
is not assumed, in which case the major symmetry ia O j/3 from ( I3.3P is lost; moreover, 
the theory could be considered more general if the explicit dependence on structure 
tensors such as r is omitted. The generality considered here, however, is sufficient for 
our purpose. 

For an incompressible material the incremental incompressibility can be written in 
either of the forms 

trfTF- 1 ) = trL = 0, divu = 0. (3.4) 
In the absence of body forces the incremental motion is governed by the equation 

DivS = p r x itt , (3.5) 

where p r is the mass density in B r and a subscript t following a comma signifies the 
material time derivative, i.e. the time derivative at fixed X, so that x t = u t is the 
particle velocity and x tt = u tt the acceleration. The incremental counterpart of the 
rotational balance equation ( 12. 4 p is 

FS + FS = S T F T + S T F T , (3.6) 

or, on use of ( 12. 6p ? and (13. ip . 

FS + JL<r = S T F T + J<tL t . (3.7) 

On introducing the 'push forward' So of S, defined by So = J -1 FS, and the corresponding 
push forward Ao of the elasticity tensor such that So = AoL (or AqL + pL — pi in the 
case of an incompressible material) we may express ( 13. 7p in the form 

AqL + L<t = (AqL) t + <tL t (3.8) 

for an unconstrained material, and as 

AoL + L(a + pl) = (AoLf + (a + pl)L T , (3.9) 

for an incompressible material. In component form A® is related to A via 

JAopiqj Fp a FqfiA a if3j , (3.10) 

with J = 1 in the incompressible case. 

For details of the background o n the theory of increm ental deformations superimposed 
on a finite deformation we refer to lOgdenl (119841 120071 ). for example. 
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From (13. 8 p we deduce that tr[(^4oL + Lcr)W] = for any skew symmetric second- 
order tensor W, and since this and the corresponding result from (I3.9P holds for arbitrary 
L we deduce further that 



A)W = -crW, AqW 



pI)W, 



(3.11) 



for an arbitrary skew symmetric W for unconstrained and i ncomp r essibl e materials, re- 
spectively. A formula equivalent to (I3.11l) i was derived by iHogerl (119861 ) by a different 
method and stated explicitly as the formula (2.2.2) in the latter paper for a residually- 
stressed but undeformed configuration. Note that the deriv ation herein does not require 
the existence of a strain-energy function. Another result in IHogerl (119861 ) is recovered by 
choosing L to be symmetric, which leads to 



;[AE-(AE) r 



1 



a-E - Ecr), 



(3.12) 



for any symmetric E. This relation applies for bo th unconstrai ned and incompressible 
materials and corresponds to the formula (2.2.3) in IHogerl ( 119861 ). the right-hand side of 
which contains a sign error. 

Note that as well as possessing the major symmetry Ao P i q j = Ao q j P i, which follows 
from (I3.3P and (13.101) . A® has the property 



for a compressible material, and 



A, 



&ij @pq i 



Aoipqj ~\~ 5ij((J 



V\ U PQ 



p5. 



pqfi 



(3.13) 



(3.14) 



for an incompressible material. These formulas follow from (I3.8P and (13.91) . respectively. 
Correspo nding formulas for a pre-stre ssed material in the a bsence of residual stress were 
given by IChadwick and Ogdenl (1197ll ) and IChadwick! (119971 ) . 

In Eulerian form, i.e. in terms of So and u, the incremental equation of motion (13. 5p 
becomes 

divS = pu tt , (3.15) 

where p = p r J~ x is the mass density in B. The (Cartesian) component form of this 
equation may be written, for an unconstrained material, as 



(AopigjUj^q) fP — pUijt, 

and for an incompressible material as 

{A piqjU jt q) t p - p ti + P tj U jti = p r Ui jtt , With U^i = 0. 



(3.16) 



(3.17) 



Suppose now that the material properties are homogeneous, implying, in particular, 
that the initial stress r is uniform, and that the deformation in B is homogeneous. Then, 
A., A$ and p are constant, and the equations of motion (I3.16P and (13.171) reduce to 



A 



Opigj^j ,pq pUi^tti 



(3.18) 
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and 

AopiqjUjpq pi = pUijt, Ui t i = 0, (3.19) 

respectively, and we note that p = p r in the latter case. 

As already indicated, we are considering A. to depend on the deformation through 
F, on the initial stress r, and on any material symmetry present in the configuration 
B r . In what follows we consider the intrinsic properties of the material to be isotropic 
relative to B r so that the only source of anisotropy is the initial stress. In order to 
account for the initial stress in the expressions for the stress and elasticity tensors we 
base the development in the following sections on scalar invariants involving the right 
Cauchy-Green deformation tensor C and the initial stress tensor r. 



4 Invariant-based formulation 

4.1 Invariants of C and r 

The strain-energy function W of an initially stressed material depends on r as well as 
the deformation gradient F, and, by objectivity, it depends on F only through C = F T F. 
For a material that is isotropic relative to B r in the absence of initial stress, W can be 
treated as a function of three independent invariants of C, which are commonly taken to 
be the principal invariants defined by 

Ji=trC, h = |[(trC) 2 -tr(C 2 )], J 3 = det C. (4.1) 

In the configuration B r these reduce to I\ = I2 = 3, J3 = 1. If there is an initial stress r 
in B r then, in general, the material response relative to B r depends also on the invariants 
of r and the combined invariants of r and C. A possible set of independent invariants 
of r, including those that depend on C, is 

trr, tr(r 2 ), tr(r 3 ), tr(rC), tr(rC 2 ), tr(r 2 C), tr(r 2 C 2 ). (4.2) 

Together, there are thus at most 10 independent invariants of C and r in general, and 
for an incompressible material we have I3 = 1 and hence there are at most 9 indepen- 
dent invariants in this case. The number of independent invariants is reduced when the 
dimension of the considered problem is reduced from three to two, as for a plane strain 
setting, for example. In the reference configuration B r the fourth and fifth invariants 
listed in (14. 2p reduce to the first, and the sixth and seventh to the s econd. For re levant 



background on invariants of tensors we refer to ISpencerl (ll97lf ) and [Zheng ( I1994I ). For 
subsequent reference we introduce the notation 

J 6 = tr(rC), J 7 = tr(rC 2 ), J 8 = tr(r 2 C), J 9 = tr(r 2 C 2 ), (4.3) 

for the invariants that depend on both C and r, noting that the notations I4 and 1$ are 
not used since these tend to be associated with a preferred direction in a transversely 
isotropic material. 

Although not immediately obvious, it is worth noting that invariants such as tr (CrCr) 
and tr (CrC 2 r) may be expressed in terms of the invariants listed above. This may be 
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shown by first applying the Cayley-Hamilton theorem to C + At for an arbitrary scalar 
A. The coefficient of A yields the identity 



C 2 r + CrC + rC 2 - h{Cr + rC) - (trr)C 2 

- ft (trr) - J 6 ]C + J 2 r - [/ 2 (trr) - IJ 6 + I 7 ]l = 0. (4.4) 

Then, by multiplying by r, for example, and then taking the trace of the result leads to 

tr(CrCr) = 2hh - 2/ 9 - 2/ 1 / 6 (trr) + 2/ 7 (trr) + l\ + J 2 [(trr) 2 - tr(r 2 )], (4.5) 

which shows how tr(CrCr) depends on the other invariants. 

The expressions for the stress and elasticity tensors given in (12. 6p . (12. 9p and (13. 3 p 
require the calculation of 



dF ^ dF 



and 



d 2 W v^ TTr d 2 Ii Tir dli dh ttns 

where we have used the shorthand notations Wi = dW/dli, Wy = d 2 W/dIidIj, i,j G X, 
and I is the index set {1, 2, 3, 6, 7, 8, 9} (or {1, 2, 6, 7, 8, 9} in the case of an incompressible 
material). Although their derivatives with respect to F vanish and do not appear in the 
above expressions, the invariants trr, tr(r 2 ) and tr(r 3 ) may nevertheless be included 
in the functional dependence of W. In the following we give explicit expressions for the 
stress tensors and the elasticity tensor based on these summations. 



4.2 Stress tensors 

For an unconstrained material the strain-energy function W depends on the seven deform- 
ation dependent invariants Ji, I2, 13, ■ ■ ■ , Ig together, in general, with trr, tr(r 2 ) and 
tr(r 3 ). From (12. 6p and ( 14. 6 p the nominal and Cauchy stresses are given by 

iex 

The corresponding expressions for an incompressible material are 

S = ^ - pF-i = £ wB - pF-\ <r — FS. (4.9) 

iex 

The required expressions for dli/dF are listed for convenience in Appendix A. In partic- 
ular, these enable the Cauchy stress to be expanded, so that 

Jo- = 2WkB + 2W 2 (h~B - B 2 ) + 2/3W3I + 2W & Y1 

+ 2ty 7 (£B + BE) + 2W 8 SB~ 1 S + 2W 9 (SB^ 1 SB + BSB ^E), (4.10) 
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where the notation X = FrF T has been introduced and we recall that B = FF T is the 
left Cauchy-Green tensor. For an incompressible material (I4.10|) is replaced by 

a- = 2WiB + 2W 2 (hB - B 2 ) - pi + 2W/ 6 £ 

+ 2W 7 (^B + BE) + 2W 8 £B- 1 £ + 2W / " 9 (SB" 1 SB + BEB^E), (4.11) 

with I 3 = 1. 

If we evaluate (14.10]) in the reference configuration, it reduces to the appropriate 
specialization of (12. 7p . namely 

t = 2(Wi + 2W 2 + W 3 )I + 2(W 6 + 2W 7 )t + 2(W 8 + 2W 9 )t 2 , (4.12) 

where all Wi, i G X, are evaluated for Ii = I 2 = 3, J 3 = 1, J 6 = J 7 = trr, 1$ = I9 = tr (r 2 ). 
This indicates that we should set 

Wi + 2W 2 + W 3 = 0, 2(W / 6 + 2W 7 ) = l, iy 8 + 2iy 9 = (4.13) 

there. In fact, if we require ( I4.12p to hold for all r then the conditions in (I4.13P necessarily 
follow. The corresponding reduction for an incompressible material yields 

t = (2Wx + 4W 2 - p (r) )I + 2(W 6 + 2W 7 )t + 2(W 8 + 2W 9 )t 2 , (4.14) 
which specializes (I2.10p . and the counterpart of ( I4.13P is 

2W X + 4W 2 - p {r) = 0, 2(W 6 + 2W 7 ) = 1, W 8 + 2W 9 = 0, (4.15) 
evaluated in the reference configuration. 

4.3 The elasticity tensor 

Next, we note that the elasticity tensor A. is given by 

d 2 W v^ TI , d 2 h TXr dh dlj , A „. 

A =imF^ w 'sfsf + £ » a#- < 4 - 16 > 

This requires expressions for the second derivatives of the invariants, which are collected 
together for reference in Appendix A in component form. Here we give the component 
form of v4o: 

JAopiqj = 2{Wi + IiW 2 )B pq 5ij + 2W 2 [2B pi B qj — B iq B jp — B pr B rq 5ij — B pq Bij] 

+ 2I 3 W 3 (25 ip 5j q — SigSjp) + 2W e T lpq 5 i j + 2W 7 {T lpq B i j + Sp,.!?^^ + B pr Yj rq bij 

+ ^ijB pq + T, p jBi q + + AWnBi p Bj q 

+ 4W 22 (IiB ip — B ir B rp )(IiB,j q — Bj S B sq ) 
+ 4/ 3 2 H / 335 i p5 J g + AWi 2 (2IiB ip B jq — B ip B jr B rq — B jq B ir B rp ) 
+ AI 3 W 13 (B ip S jq + B jq 5 ip ) + 4/3^23 [^(-Bip^g + B jq 5 ip ) — 5 ip B jr B rq — 5 jq B ir B rp ] 
BjqT.ip) + 4:Wi 7 [B ip (J2j r B rq + Bj r T, rq ) + Bj q {T, ir B rp + B ir Y> rp )\ 
+ W 26 [(hB ip -B ir B 

rpj^jq + \I\Bjq Bj r B rq )T lip ] 

+ 4W 27 [(IiB ip — B ir B rp )(Ej S B sq + Bjs'Zsg) + (hBj q — Bj r B rq )(E is B sp + B is E sp )} 
+ 4I 3 W 36 {5 ip Z jq + 5 jq E ip ) + AI 3 W 37 [8 ip (Ej r B rq + Bj r Yj rq ) + 5j q (Yli r B rp + 5j r S r p)] 

+ 4H / 66Si p Sjg + 41^67 [Sjp(S Jr i? rg + Bj r T, rq ) + YtjqiYjirBrp + I?j r E r p)] 

+ 4W 77 (T, ir B rp + Ej r .E r .p)(S JS S S9 + Bj S E sq ), (4.17) 
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and, in order to save space, terms involving derivatives with respect to I 8 or I 9 have not 
been included. Note, however, that it is straightforward to include these terms, just by 
repeating the terms involving derivatives with respect to Iq and I 7 , replacing the indices 
6 and 7 by 8 and 9, respectively, and replacing S = FrF T by Fr 2 F T . When there is 
no initial stress the above formula reduces to the resul t for a prestressed i sotro pic elastic 
solid, as given in, for e xamp le, t he classica l pape r by lHaves and Rivlinl ( 1l96ll ); see also 



Toupin and Bernstein! (119611 ) and llruesdelll (119611 ) for related basic theory. 



The corresponding formula for an incompressible is obtained by setting J = 1 and 
J 3 = 1 and omitting the terms involving derivatives of W with respect to J3, and yields 



Biq Bjp 



Bp r B r q5ij 



Bp r ^j r q5{j 



^ij Bpq 



BpqBij] 
^pjB'iq 



^iqBpj] 



B j sB S q) 



A 0piqj = 2(Wt + hW 2 )Bp q 5 lj + 2W 2 [2B pi B qj 
+ 2WoLp q 5ij + 2W-j\L pq Bij + I]p r B r q5i 
+ 4W n B ip B jq + m 22 [I 1 B i p - B ir B rp ){hB jq 
+ 4W 12 (2I 1 B ip Bj q — BipB jr B rq — B jq B ir B rp ) 

+ 4Wi 6 (-BipEjg + BjgTjip) + AWnlBip^jrBrg + BjrT^rq) + Bj q (E ir B rp + B ir T, rp )\ 

+ ^W 2& [(IiBip — B ir B rp )E jq + (hB jq — B jr B rq )E ip \ 

+ 4W 27 [(IiB ip — B iT B r p){Yjj S B sq + Bj S T, sq ) + {I\Bj q — Bj r B rq ){Tj is B S p + S is E sp )] 

ip^jq + 4W / 67pip(E J >_B r(J + BjrTjrq) + Yjj q {Tji r B rp + Bi r Ti rp )] 

+ 4W 77 (E ir B rp + B ir E r p)(Ej S B S q + 5 JS S sg ). (4-18) 

When evaluated in the reference configuration, the moduli (I4.17p . with the missing 
terms included, reduce to 

•A-Opiqj ^1 ifiij^pq ^iq^jp) ^2&ip3jq 8%jT pq 01\^ijr pq ~\~ Sp q T{j -\- S{ q Tjp -\- SjpTi q ^j 

+ h^ipTjq + 5 jq Tip) + (3 3 T ip T jq + ^(SijTpkTkq + S pq T ik T k j + Si q T jk T kp + 5j p T ik T kq ) 
+ l2{Si p T jk T kq + 5j q T ik T kp ) + ^(TipTjkTkq + T jq T ik T kp ) + ^Aji k T k pTj k T kq , (4.19) 

where the a's, /3's, and 7's are defined by 

olx = 2(Wi + W 3 ), a 2 = 4(Wn + AW 12 + AW 22 + 2W 13 + AW 23 + W 33 ) - 2a u 

= 2W 7 , p 2 = 4(Wie + 2W^i 7 + 2W 26 + 4^ 27 + W 3e + 2W 37 ), 
03 = 4(^66 + 4^ 6 7 + 4^7), 71 = 2W 9 , 

72 = 4(W 18 + 2W 19 + 2W 28 + 4^29 + ^38 + 2W 39 ), 

73 = 4(W m + 2W m + 2^8 + 4^79), 74 = 4(jy 88 + AW 89 + 4W 99 ), (4.20) 



all derivatives being evaluated in the reference configuration and use having been made 
of the connections (14.13)) . Note that in general the expressions (14.20)) may depend on the 
invariants trr, tr(r 2 ) and tr(r 3 ). 

It is worth noting here that when referred to axes that coincide with the principal 
axes of r, the only non-zero components of (14. 19ft are given by 



Aoiiu = 2a x + a 2 + (1 + 4ft + 2(3 2 ) Tl + /3 3 r 2 + 2(2 7l + 72 )r 2 + 2 73 rf + 7 4 r, 
Aoujj = a 2 + p 2 (Ti + Tj) + /3 3 Tii~j + liijl + rf) + + r^Tj + 74^?, 
Aoijij = a 1 + r i + Pi{Ti +Tj) + 7i(^ 2 + rf ) = A 0ijji + r<, 



1 

i 1 



(4.21) 
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where there are no sums on repeated indices, % ^ j, and n (i — 1, 2, 3) are the principal 
values of r (in general, 15 non-zero components of Ao in total). It is interesting, but 
not surprising, to note that for a prestressed isotropic elastic material the only non-zero 
comp onents of Aq are also Aouu, Aoujj, Aoijij and Aoijji, i ^ j (see, for example, lOgdenl . 



198J). 

For an incompressible material there is a further reduction of (14.191) . leading to 

Aopiqj OLy Spg -\- &iq$jp) ~\~ &ij7~pq fil(3ijT~pq ~\~ $pqT~ij ~\~ &iqT~jp ^jpTiq} 

+ ^{^ipTjq + SjqT ip ) + faT ip T jq + ^l(5ijT pk T kq + S pq T ik T kj + S iq T jk T kp + Sj p T ik T kq ) 
+ l2{5ipTj k T kq + S jq T ik T kp ) + ^{T ip T jk T kq + T jq T ik T kp ) + r )4 : Ti k T kp Tj k T kq , (4.22) 

since the term involving a 2 may be dropped because it vanishes from the expression for 
the incremental stress by virtue of the incompressibility condition. Here ot\, fli, (3 3 , 71, 73 
and 74 are unchanged, while $2 reduces to 4(Wi6 + 2Wn + 2H^6 + 4H^ 2 7) and 72 reduces 
to 4(Wis + 214^19 + 214^28 + 41^29)- The terms foSjqTip and ^SjqTikTkp also vanish from the 
incremental stress for the same reason but are retained here since otherwise the major 
symmetry would be lost. Note that the terms involving 5j P in (I4.22p do not contribute 
to the equation of motion, again by virtue of the incompressibility condition. 

5 The effect of initial stress on infinitesimal wave 
propagation 

Now consider incremental motions in an infinite homogeneous medium subject to homo- 
geneous deformation and/or homogeneous initial stress. From Section [3] we recall that 
the equation of incremental motion is 

•Aopiqj u j,pq = P u i,tt? (5-1) 

for a compressible material, and 

AopiqjUjpq pi = pUi^tti Ui^i = 0, (5-2) 

for an incompressible material. 

Consider a homogeneous plane wave of the form 

u = m/(n • x — vt), (5.3) 

where m is a fixed unit vector (the polarization vector), / is an arbitrary function of the 
argument n • x — vt, n is a unit vector in the direction of propagation, and v is the wave 
speed. For an incompressible material we also set 

p — g(n ■ x — vt), (5.4) 

where g is a function related to /. Substitution into the equation of motion (15. ip (after 
dropping /", which is assumed to be non-zero) leads to 

AopigjUprigrrij = pv 2 mi, (5.5) 
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for a compressible material, and into ( 15. 2 p 

Aopiqjnprigmjf" - g'm = pv 2 rriif", m^i = 0, (5.6) 

for an incompressible material, where a prime on / or g indicates the derivative with 
respect to its argument. By multiplying the first equation in (15. 6p by m and using the 
second equation, we obtain 

g' = AopigjUpngmjUif", (5.7) 
and upon substitution back in (I5.6P for g 1 and again dropping /" we arrive at 

(A Q pi qj npn q - AQ pkq jn p n q n k ni)mj = pv 2 mi. (5.8) 

It is convenient now to use the acoustic tensor Q(n), whose components are defined 

by 

Qijip) = AopigjUpUq. (5.9) 
Then equations (15.51) and (15. 8p can be written in the compact forms 

Q(n)m = pv 2 m, (5.10) 

Q(n)m = pi; 2 m, m ■ n = 0, (5-H) 

for unconstrained and incompressible materials, respectively, where, noting that Im = m, 
with I = I — n <g> n defined as the projection operator on to the plane normal to n, the 
symmetric tensor Q(n) is defined as 

Q(n) = IQ(n)I. (5.12) 

This is the projection of Q(n) on to the plane normal to n. The symmetrization of the 
acous t ic ten sor in the incompressible c a se is attrib uted to Hayes in the Ph.D. thesis of 
Scottl (11974^ : see also IScott and Haves! (119851 ) and iBoulanger and Havesl (119931 ) . 



For any given direction of propagation n we now have a symmetric algebraic eigenvalue 
problem for determining pv 2 and m, which is three-dimensional in the compressible case 
and two-dimensional in the incompressible case. Because of the symmetry there are 
three mutually orthogonal eigenvectors m corresponding to the direction of displacement 
in the compressible case and two in the incompressible case, corresponding to transverse 
waves. The values of pv 2 , three and two, respectively, are obtained from the characteristic 
equations 

det[Q(n) - pv 2 I] = 0, det[Q(n) - pv 2 I] = 0, (5.13) 

for compressible and incompressible materials, respectively. 
If m is known then pv 2 is given by 

pv 2 = [Q(n)m] • m, (5.14) 

which applies for either compressible or incompressible materials. 



■""We are grateful to a reviewer for bringing this reference to our attention. 
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To examine plane waves in detail we need to calculate Q(n) for the forms of «4o given 
in Section fl~3l We now do this explicitly for the reference configuration B r . From (I4.19P 
and (I5.9p we obtain 

Q(n) = [«! + (1 + /3i)(n ■ rn) + 7l (n ■ r 2 n)] I + far + 7iT 2 

+ («i + a 2 )n ® n + + (3 2 )(n <g> rn + rn <g> n) + /3 3 rn (g) rn 

+ (71 + 72) (n <g> T 2 n + r 2 n <g> n) + 73(711 g) r 2 n + r 2 n g) rn) + 74 r 2 n <g> r 2 n 

(5.15) 

for compressible materials, and from (I4.22|) . ( 15. 9p . and (15.121) . we obtain 

Q(n) = [ ai + (1 + /3i)(n ■ rn) + 7l (n ■ r 2 n)] I + ftlrl + 7 iIr 2 I 

+ /3 3 Irn <g> Irn + 73 [Im <g> Ir 2 n + Ir 2 n <g> Im] + 74 lT 2 n <g> Ir 2 n (5.16) 

for incompressible materials. Note, in particular, the explicit nonlinear dependence of 
the acoustic tensors on the initial stress r (not forgetting that the a's, /3's, and 7's may 
themselves be nonlinear in r). 

Example 1: Longitudinal and transverse waves in compressible materials We 

consider the possibility of a pure longitudinal wave in the reference configuration (initial 
stress, no pre-strain), so that m = n. From (15.101) and (I5.15P we obtain 

(a — pv 2 )n + bm + CT 2 n = 0, (5-17) 

where the constants a, b, c are given by 

a = 2e*i + a 2 + (1 + 2fr + f3 2 )(n ■ rn) + ( 7i + 7 2 )(n ■ r 2 n), 
b = 2ft + 02 + /? 3 (n • rn) + 73 (n • r 2 n), 

c = 2 7 i + 72 + 7 3 (n ■ rn) + 7 4 (n ■ T 2 n). (5.18) 

Several distinct cases arise for which a longitudinal wave can propagate provided the ma- 
terial parameters and the initial stresses are such that pv 2 > 0. First, for the degenerate 
case in which b = c = we have pv 2 = a and a longitudinal wave can propagate in any 
direction. Second, if c = and b ^ then n is an eigenvector of r. Recalling that 
Tj, 2 = 1,2, 3, are the principal values of r, a longitudinal wave can then propagate along 
the z-th principal direction, % = 1,2,3, with speed v given by pv 2 = a + bTi. Third, if 
c^O then we refer (I5.17P to the principal axes of r, leading to 

{a- pv 2 + br l + cr 2 )n l = 0, i = l,2,3. (5.19) 

There are now different possibilities depending on the multiplicity of the eigenvalues of r. 
If Tj = r say, i — 1, 2, 3, then pv 2 = a + br + cr 2 and a longitudinal wave can propagate in 
any direction. If two of the eigenvalues are equal but distinct from the third, say n = r 2 = 
t ^ T3, then a longitudinal wave can propagate along the 3 direction with speed given by 
pv 2 = a + br 3 + ct 2 or in the principal (1,2) plane with speed given by pv 2 = a + br + cr 2 . 
There is also a special case in which if b + c(r + r 3 ) =0 then a longitudinal wave can 
propagate in an arbitrary direction with speed pv 2 = cl—ctt^. Finally, if the eigenvalues Tj 
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are distinct then, for propagation in the i-ih principal direction, the wave speed, expanded 
in full, is given by 



pv 2 = 2at + a 2 + (1 + 4ft + 2ft)r, + (ft + 3 7l + 2 72 )r 2 + 2 73 r? + 74 r/. (5.20) 

It is also easy to show that propagation is possible in any direction within an 
principal plane if b + c[ji + Tj) = 0, in which case the formula for the wave speed is 

pv 2 = a — CTiTj. 

Of course, a longitudinal wave is in general accompanied by a pair of transverse waves 
with mutually orthogonal polarizations. In the case where a pure longitudinal wave exists 
and n is aligned with the principal axis of r, say, with i — 1, 2 or 3, it follows that the 
expression (15 . 1 5[) for the acoustical tensor reduces to 

Q(eO = [a x + (1 + ft)ri + 7 irf]I + ftr + 7l r 2 

+ [ ai + ct 2 + 2(ft + /3 2 )n + (ft + 2 7l + 2 72 )r 2 + 2 73 rf + l4 rf\ ® e,. (5.21) 

Hence, for a corresponding transverse wave with polarization m and speed v, we have, 
on use of (I5.10p . 

[a x + (1 + Pi)n + 7i7f ]m + Arm + 7l r 2 m = pw 2 m. (5.22) 

In general it does not follow that m is also a principal axis of r. However, in the special 
case in which m is also a principal axis of r (in the plane normal to ej), let this correspond 
to principal initial stress Tj, j ^ i. Then, the wave speed, denoted v^, associated with 
m = ej is given by 

pv f. = Ti + ai + ft( T . + Tj ) + 7l ( r 2 + r 2). (5.23) 

The difference 

P(vl - v%) = r % - r v (5.24) 



recalls the universal relationship est ablished by iMan and Lul (119871 ) for the general lin 



ear theory. Note, however, that in Man and Lul (119871 ) the material was taken to be 



orthotropic with the principal axes of r coinciding with the axes of orthotropy, whereas 
here the orthotropy is associated with r itself and the underlying material is isotropic. 
The principal axes are eigenvectors of Q in each case although the values of the compo- 
nents Qij differ. 

The situation is somewhat different if one starts by considering the existence of a 
pure transverse wave. There is then no general result showing that such a wave must be 
travelling and/or polarized along a princip al axis of initial stres s. We do not include here 



the details of this case, and refer instead to IMan and Lul ( 119871 ) for an example of a wave 



polarized along one principal direction of stress and propagating in any direction in the 
plane normal to the polarization. 

In all the above examples, for a real wave to exist it is necessary that pv 2 > 0, which 
puts restrictions on the parameters and initial stress components involved. In general 
this requirement may be stated in terms of the strong ellipticity condition, i.e. 

m ■ [Q(n)m] > for all non-zero vectors m, n. (5.25) 

This also applies in the incompressible case, subject to the restriction m • n = 0. 



15 



Example 2: Connections with the results of Man and Lu Suppose now that in 
(I5.15P the coefficients a%, 0:2, ft, ft are constants and only the terms that are linear in the 
initial stress r are retained. In particular, the terms with coefficients ft, 7$ (i = 1, 2, 3, 4) 
are omitted. Then, referred to the principal axes of r, the components of Q(n) can be 
written compactly as 



Qu = cm + (ai + a 2 )n 2 i + (1 + ft)(n ■ rn) + ftr* + 2(ft + ft)T;n; 



1,2,3, 



Q 



["i + «2 + (ft + ft) fa + r^Kn,,-, z 7^ j, 



(5.26) 



where r, (« = 1,2,3) are the principal values of r. These formulas are consistent with 
correspondi ng expressions for the components of the acoustic tensor (or Christoffel tensor) 
obtained by Man and Lu f 19871 ) [their eq. (18)] although this is not immediately clear 
since the notatio n differs. The corr espondence can be established by noting that the 
tensor C used by lMan and Lul (119871 ) has components that are connected to those of *4o 
by Cijki = Aoijki — bjiTik- 



E xample 3: Connect ions with the classical theory of Biot In the classical theory 
of iBiotl (119391 Il940l . 119651 ) the incremental stress may be written (in the present notation) 



as BL or in index notation as Bpi q jUj >q , where the coefficients B P i q j satisfy the relations 



B = B — B 



B — B 



(5.27) 



the latter following from the existence of a strain-energy function. It was shown in 
Ogden and Singhl (1201 if ) that the general connection between Ao P iqj and B P i q j may be 
written in the form 



A 



Opiqj 



B 



PK-U 



IS t ■ - -S r ■ - -S -t ■ + -S- t + 5 -T ■ 

2 PJ 1 % VQ 13 9 Q1 1 V3 ~ 9 U 13 ' VQ ~ U Q3 1 PI 



2 U P1 ' V 



2 w Q l ' P0 



2 Wl 3 ' PI 



>qj 'pi- 



(5.21 



It is worth noting that this connection does not depend on the existence of a strain-energy 
function and is independent of any material symmetry. For the general expression (I4.19P 
to reduce to the Biot form for isotropic response the material parameters in (I4.19P must 
be specialized to 



«i = /i, 



«2 



A. 



ft 



'1/2, ft 



1. 



(5.29) 



where fi and A are the Lame moduli, and the terms which are of order higher than 1 in 
t must be neglected. This corresponds to B P i q j having the form 



B 



H(Sij8 pq + S qi S pj ) + \5pi5qj + S pi T ( 



<IJ- 



(5.30) 



which satisfies the conditions (15.271) . 

Clearly, the formula ( 15.301) relies on the special identifications ( I5.29p . and is very much 
a specialization of the more general theory discussed here. In effect, if one was to linearize 
the general expression of the elastic moduli (I4.19P with respect to the initial stress, then 
the expansions 



Cti = /i + &iT, Oi2 = A + &2T, 



ft =01, ft = ft 



(5.31; 
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would be required, where "r" is a term of first order in r, and the scalars with a hat are 
constants, which cannot be determined a priori and must be measured experimentally. 
In conclusion, six material constants (A, p, a±, a 2 , fii, (3 2 ) are required to describe the 
isotropic response of a compressible material with a sma ll initial stress, not just the two 
Lame coefficients as Biot implied. We refer to iMan fll998h for a discussion of this point in 



relation to Hartig's law and a proof that four additional constants are needed to describe 
the elastic response of currently known real materials. 

Example 4: Second-order elasticity Here we consider an isotropic elastic solid 
without initial stress and its specialization to second order in the components of the 
Green strain tensor E = |(C — I), where C is again the right Cauchy-Green deformation 
tensor. An appropriate fo rm of the e l astic strain-energy function, which is correct to 



second order, is that due to lMurnaghanl (119371 ). When expressed in terms of the invariants 



of C given by (14. ip the Murnaghan energy function has the form 
W = ^(ii - 3) 2 + ^(I 2 - 2ii - 2I 2 + 3) 

+ ^(h ~ 3) 3 + -(/i - 3)(/ 1 2 - 3/ 2 ) + -(h -I 2 + h- 1), (5.32) 

where A and p are the classical Lame constants of linear elasticity and I, m, n are the 
second-order constants of Murnaghan. For this energy function we have W 22 = W13 = 
1^23 = ^33 = and since we are not including initial stress the expression (I4.17P reduces 
to 

JAopiqj = 2{W\ + I\W 2 )B pq 5ij + 2W 2 [2B P iB q j — B iq Bj P — B pr B rq 5ij — B pq Bij] 
+ 2hWz(25 ip 5 jq - 5 iq 6 jp ) + 4W u B tp B jq 

+ AW 12 (2I 1 B ip B jq — B ip B jr B rq — B jq B ir B rp ), (5.33) 

and the remaining coefficients Wi, W 2 , W 3 , Wu and W± 2 are simply obtained from f !5.32[) . 
Then, with this specialization, the wave speed v is given by 

pv 2 = AopigjUpngmiiTii. (5.34) 

In working with second-order elasticity the corrections to the classical longitudinal and 
transverse wave speeds are obtained to first order in E. We have C = I + 2E and 
l\ = 3 + 2E, exactly, which we use together with the linear approximations I 2 ~ 3 + 
I 3 ~ 1 + 2E, B ~ C = I + 2E, and we also note that p ~ p r (l — E), where E = trE. 
To the first order in E we then obtain 

Wi = p + |n + |(A + 2 fx + 2m) E, W 2 = -\p - \n - \mE , W 3 = \n, 
W 11 = \{\ + 2p + Am) + \{l + 2m)E, W 12 = -\m. (5.35) 

After some manipulations the above approximations enable us to obtain, to the first 
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order in E, 



J <Aopiqj ft\3ijfipq ^iq^jp) ~\~ X5{p5jq ~\~ ^^i^^ijEpq + 5pqE{j -\- 5{qEjp + 5jpE^q) 
+ ^{E5ij5 pq + 2S ip Ejq + 25j q E ip ) + 2lES ip 5jq 

+ m[E(5ijSpq + 5 iq Sjp - 25ip5 jq ) + 2(6 ip E jq + S jg E ip )] 

-\- -^Tl\&ijEp q -\- 5 p qEij -\- SiqEjp -\- SjpE^q 2SipEjq 25jqEip 

- E(5ij5 P q + 5iq5 jp - 25i P 5jq)], (5.36) 

and hence, for any m, n pair satisfying (j5.10|) . the wave speed is obtained via 

p r v 2 = fi + (fi + A)(m • n) 2 + 2/i[2(n • En) + (m • Em) + 2(m • n)(m • En)] 
+ X[E + 4(m • n)(m • En)] + 2lE(m ■ n) 2 
+ m{E[l - (m ■ n) 2 ] + 4(m • n)(m • En)} 

+ |n{(n ■ En) + (m ■ Em) - 2(m ■ n)(m ■ En) - E[l - (m ■ n) 2 ]}. (5.37) 

For a longitudinal wave with m = n this reduces to 

p rV 2 = A + 2/i + (A + 21) E + (10// + 4A + 4m) (n ■ En), (5.38) 

and for a transverse wave with m • n = and m x n = 1, 

p r . v 2 = fj, + 2/i(n • En) + (A + 2/x + m)E - |(4/i + n)(l ■ El), (5.39) 

where we have made use of the connection 1 ■ El + m ■ Em + n ■ En = E. When n is 
specialized to the axis e% and m is taken to be e 2 and e 3 in turn, wher e gj , e? and e 3 are 



princi pal axes of strain, then these results agree with those obtained by lHughes and Kelly 



(jl953|) [their eq. (11)]. 

It is straightforward to show that a longitudinal wave actually exists if either n is 
an eigenvector of E or the elastic constants satisfy A + 2[i + m = 0. For the existence 
of a transverse wave there are more options. First, in the special case when the elastic 
constants are such that A + 2/i + m = and 4/x + n = 0, a transverse wave exists for 
any direction of propagation n. Second, if 4(A + /i + m)— n = and 4/i + n ^ then 
a transverse wave exists only if m is an eigenvector of E. Third, if 4/i + n = and 
A + 2/i + m ^ then a transverse wave exists for any n for which Em lies in the plane 
normal to n. Finally, if none of these very special cases apply then the existence of a 
transverse wave requires that m be an eigenvector of E. 



Note that the case of incompressible materials is treated elsewhere; see lDestrade et al. 



(120101 ) and references therein. 



Example 5: A simple nonlinear model The above examples have involved different 
specializations that consider the strain and/or the initial stress to be small. In this final 
example we consider a finite deformation from a configuration that is subject to an 
initial stress that is nonlinear. For this purpose it suffices to adopt a simple prototype 
incompressible form of strain-energy function, which is taken to be 

W = \n{h - 3) + 1/2(J 6 - trr) 2 + |(J 6 - trr), (5.40) 
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where the first term represents the classical (isotropic) strain-energy function of a neo- 
Hookean material with shear modulus /x, while the terms in Iq capture the effect of initial 
stress and satisfy (14. 15|) q . The constant Jx has dimensions of [stress] -1 . We recall that 
Iq = tr(Cr), which may also be written as trS. 

From (14.1 II) the Cauchy stress is then obtained as 

CT = /jB-pI+[l + 2/i(I 6 -trT)]E, (5.41) 

and from (14. 18ft the components of »4o reduce to 

Myiqj = I^BpgSij + [1 + 2//(/ 6 - trr)]E M <% + AfxT, pi E qj . (5.42) 

The corresponding components of the acoustic tensor are obtained from (15.91) . and it is 
convenient to write it in tensor notation as 

Q(n) = al + /3Sn <g> En, (5.43) 

where for compactness of representation we have introduced the notations 

a = /i(n-Bn) + [l + 2/i(/ 6 -trT)](n-Sn), ,9 = 4/2. (5.44) 

Since we are considering an incompressible material we form the projection of Q(n) 
on to the plane normal to n according to (I5.12p . This is 

Q(n) = al + /3ISn <g> ISn. (5.45) 

Next we solve A5.13[) q and find that the two solutions for pv 2 are 

pv 2 = a, pv 2 = a + /3(n- Em) 2 , (5.46) 

and equation (15. lip becomes 

am + /3(n • Sm)ISn = pv 2 m, m ■ n = 0. (5-47) 

Assuming that (3 ^ 0, for the solution pv 2 = a this requires that either ISn = or 
n • Sm = 0. If ISn = then n is an eigenvector of S and both transverse waves 
have the same speed v, given by pv 2 = a. If, on the other hand, n • Sm = and 
ISn ^ then the wave speed corresponding to polarization m is given by pv 2 = a, 
while for the second transverse wave, with polarization m' say, the wave speed is given 
by pv 2 = a + /3(n • Sm') 2 . 

As a simple illustration of the above results we now consider the initial stress to be 
uniaxial, of the form r = t\B\ <E> ei, where ei is a principal axis of B corresponding to 
principal stretch Ai. Then, S = TiAfei ® ei. It follows that n • Sm = and 

pv 2 = ( fi + Tl )A 2 + 2ht 2 X 2 1 (X 2 1 - 1). (5.48) 

It now suffices to measure the wave speed for three different values of Ai to evaluate 
the material constants /t and /2, and the initial stress T\. If the deformation is achieved 
by applying a uniaxial stress along ei then the associated principal Cauchy stress o\ is 
obtained from (I5.4ip as 

ox = (p + tOA 2 -p + 2/2r 2 A?(A 2 - 1). (5.49) 

By symmetry, incompressibility and vanishing of the lateral stress it also follows from 
(I5.4ip that p = fxXi 1 . If this is used in the above then the expression for the wave speed 
can be written simply as pv 2 = a\ + /iA^ 1 . 
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Appendix A 

Al. First derivatives of the invariants 

Here we list the first derivatives of the invariants I\, I2, 13, Iq, ■ ■ ■ , Ig with respect to the 
deformation gradient. These are 

§ = 2FT ' § = 2 ^ FT " FTFFT )< § = 2/ 3F- 1 , 

^ = 2rF T , ^ = 2rF T FF T + 2F T FrF T , 

oF 3F 

91s _ Ola 



2r 2 F T , = 2t 2 F t FF t + 2F T Fr 2 F T 



from which we obtain the expressions 

F i- 2B > F§ = 2(/ lB -B^ F§ = 2/ 3 I, 

01 01 
F^ = 2S, F^ = 2£B + 2B£, 

31 31 
F-^ = 2EB- 1 £, F— - = 2SB- 1 SB + 2BSB- 1 S, 
3F 3F 

required in the expansion of Cauchy stress, where S = FrF T and B = FF T . In the 
reference configuration B r these reduce to 

2^ = 21 ^ = 41 ^ = 21 ^ = 2r ^ = 4r ^ = 2r 2 ^ = 4r 2 
3F ' 3F ' 3F ' «9F ' 3F ' <9F ' <9F 

A2. Second derivatives of the invariants 

Here we present the second derivatives of the invariants in index notation, omitting the 
details of the calculations. In the form required for the calculation of the components of 
Aq we have 

P P d 2/l _ OR X 

m qP 3F ia 3F 3P pq iv 
3 2 I 2 

FpaFqfi — <fw7r ^Il-BpqSij ~\~ ABp^Bqj IB^qB jp 1Bp r B r q5ij IBpqB^j , 

3 2 h 

F pa F ql3 — -Qjjr - 4I 3 5ipdj q - 2I 3 5 iq 5 jp , 
pa t 

fpa-Tq/3 Qp 3F j3 ~ Pq J '' 
(9 2 / 7 

FpaFqfi 2,^jpqB{j -\- 2,^2 pr B r q6ij ~\~ 2,Bp r ^j r q5{j ~\- 2,Yj{j Bpq -\- 2,^jpjB{q -\- 2,Yjq{Bj p . 

OtiaOtjfi 
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The second derivatives of Is and ig can be deduced from those of Iq and I7, respectively, 
by replacing r by r 2 , or X by SB _1 5]. In the reference configuration B r these reduce to 

d 2 I\ d 2 l2 



8F ip dF n pqi3 ' OF ip dF jq ~ pq iJ [ ip 3q iqjp ' 

8F ip dF n ~ ° lv ° 3q ° iq ° 3Jn dF ip dF jq ~ lTpqQih 
d 2 I 7 

—————— — QTp q 5ij + 2Tij5 pq + 2Tj p 5i q + 2Ti q 0j p . 

(J"i p Orjq 
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